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An effective field theory is derived for the ferromagnetic transition of diffusive electrons at 
T = 0. The static disorder which leads to diffusive electron dynamics induces an effective long- 
range interaction between the spins of the form l/r 2d ~ 2 . This leads to unusual scaling behavior 
at the quantum critical point, which is determined exactly. The crossover from this quantum fixed 
point to the classical Heisenberg fixed point should be observable in ferromagnetic materials with 
low Curie temperatures. 

PACS numbers: 64.60. Ak , 75.10.Jm , 75.40.Cx 

Quantum phase transitions have attracted substantial interest in recent years. They occur at zero temperature 
(T = 0) as a function of some non-thermal control parameter, and the relevant fluctuations are of quantum rather 
than of thermal nature. If the quantum phase transition has a classical analogue at finite temperature, then the former 
tends to be simpler than the latter in the physical dimension d = 3. For instance, the ferromagnetic-to-paramagnetic 
transition of itinerant electrons at T = as a function of the exchange interaction strength r t , has a mean- field 
like critical behavior in d = 3 Q|, since the coupling of statics and dynamics that is inherent to quantum statistics 
problems effectively increases the dimensionality of the system from d to d + z with z the dynamical critical exponent. 
In the case of clean quantum ferromagnets, z = 3, which reduces the upper critical dimension from d+ = 4 in the 
classical case to d+ = 1 M. 

If one adds quenched disorder in the form of static, nonmagnetic scatterers to the problem, then the resulting 
diffusive dynamics of the electrons changes the dynamical exponent z. Hertz has proposed that z = 4 in the 
disordered case, which would imply that <i+ = 0, that the critical behavior is otherwise unchanged, and that in 
particular the correlation length exponent v — 1/2 for all d > 0. It is now known that this cannot be true for d < 4, 
since v = 1/2 would violate the Harris criterion, which states that any fixed point (FP) with v < 2/d is unstable 
with respect to quenched disorder J^,§]. In this Letter we show that Hertz's proposition is true only for d > 6, while 
for 2 < d < 6 the disordered itinerant quantum ferromagnet shows nontrivial critical behavior which we determine 
exactly. The critical behavior is dominated by the strong static correlations that are induced by the disorder, not by 
the quantum fluctuations. Indeed, we find z = d for 2 < d < 4, which leads to d+ = 2, so that we work above the 
upper critical dimension. However, the disorder leads to additional soft modes that lead to an effective long-range 
interaction between the spins, and that are responsible for the nontrivial critical behavior. 

Our results can be summarized as follows. For the correlation length exponent v, the order parameter susceptibility 
exponent ry, and the dynamical exponent z, we find for 2 < d < 4, 

v = l/(d-2) , r) = 4-d , z = d , (1) 

and v = 1/2, r\ — 0, z = 4 for d > 4. These exponents 'lock into' mean-field values at d = 4, but have nontrivial 
values for d < 4. The magnetization, m, at T = in a magnetic field H is given by the equation of state 

tm + v m d/2 +urn 3 = H , (2a) 

where t is the dimensionless distance from the critical point, and u and v are finite numbers. From Eq. (|2a| ) one 
obtains the critical exponents (3 and S, defined by m ~ and m ~ H 1 / 6 , respectively, at T = 0, as 

= 2/(d-2) , 6 = d/2 , (2b) 

for 2 < d < 6, and (3 — 1/2, 6 — 3 for d > 6. These exponents 'lock into' their mean-field values only at d = 6. 
Equation ( paj ) applies strictly to T — 0, so t can be changed by changing the disorder or T t , but not the temperature. 
At finite temperature, we find that m obeys a homogeneity law, 

m(t, T, H) = b- f3/v m{tb 1 / v ,Tb' i ' /v , Hb^l v ) , (3) 
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where b is an arbitrary scale factor. <p — 2v for all d > 2 is a crossover exponent that describes the crossover from 
the quantum FP to the classical Heisenberg FP. The behavior of the magnetic susceptibility, Xm, follows from Eq. 
(||) by differentiation with respect to H. In particular, the exponent 7, defined by Xm ~ t ~ 7 at T = H = 0, is 
7 = (3(5 — 1) = 1 for all ci > 2. Notice that the temperature dependence of the magnetization is not given by the 
dynamical exponent, z controls the temperature dependence of the specific heat coefficient, 7y = cy/T, however, 
which has a scale dimension of zero for all d, and logarithmic corrections to scaling for all d < 4, 

lv {t,T,H) = e('i-d)\nb + lv (tb 1 ^,Tb z 1 Hb 5 ' 3 ^) . (4) 

Equations (|l|) - (^) represent the exact critical behavior of disordered itinerant quantum Heisenberg ferromagnets in 
all dimensions d > 2 with the exception of d = 4. There are several remarkable aspects of these results. First, we are 
able to obtain the critical behavior in d = 3 exactly, yet it is not mean-field like. The exactness is due to the fact that 
we work above the upper critical dimension d+ = 2. The nontrivial exponents are due to the presence of noncritical 
soft modes, or spectator modes, which lead, e.g., to the unusual term ~ v in Eq. (|2a|) . These points will become clearer 
in connection with Eqs. (Q), and (|J) below. Another point of interest is the appearance of a logarithmic correction 
to scaling in Eq. (Q). In a renormalization group context, logarithmic anomalies can arise either due to marginal 
operators, or if a set of scale dimensions fulfills some resonance condition In the present case the logarithm can 
be traced to the condition z = d, which is fulfilled for all 2 < d < 4. The behavior in d — 4 is more complicated. ||. 

The above values for the critical exponents j3 and 8 are very different from those for classical Heisenberg magnets 
H . Consequently, we expect a pronounced crossover from quantum to classical critical behavior as the ferromagnetic 
transition is approached at low temperatures. Our results can be experimentally checked in alloys such as Ni-Al or 
Ni-Ga 0. Especially the small value of 8 = 3/2 in d = 3 should be a prominent feature that could be observed for 
a fixed alloy composition as a function of magnetic field. A more thorough discussion of experimental aspects of our 
theory, and a complete account of technical details, will be given in Ref. ||. 

In the remainder of this paper we sketch the derivation of these results. Let us consider a d-dimcnsional system 
(d > 2) of interacting electrons in the presence of nonmagnetic quenched disorder. We keep the disorder strength 
sufficiently small for the system to be metallic || . With increasing strength of the exchange interaction this system 
undergoes a phase transition from a paramagnetic metal to a ferromagnetic metal. At T — this is the quantum phase 
transition we are concerned with. At any T > the asymptotic critical behavior is governed by a classical FP, but at 
low temperatures there still exists a sizeable quantum critical region that is controlled by the zero temperature FP. We 
proceed to derive an order parameter description of the quantum phase transition, along the lines of Ref. |Q. To this 
end, it is convenient to split off from the action the particle-hole spin-triplet channel, whose coupling constant we denote 
by Pj. This four-fermion interaction term we decouple by means of a Hubbard-Stratonovich (H-S) transformation 
with a H-S field M(q). All degrees of freedom other than M are then integrated out. This procedure in particular 
integrates out diffusive modes in the spin-triplet channel. These are the soft 'spectator' modes, which have important 
consequences as we will see. We obtain the partition function Z in the form Z = e~ F °l T J D[M] exp[— $[M]1 , with 
Fq the noncritical part of the free energy. A formal expansion of the Landau- Ginzburg- Wilson (LGW) functional $ 
in powers of M takes the form JjJ , 

$[M]= / M(g)[l-r tX (2) (g)]M(-g) + -L f x < 4 ) g 2 , 93 ) (M( gi ) • Mfe)) 

x (M(q 3 ) ■ M(- qi -q 2 - g 3 )) - W H M z (q - 0) . (5) 



Here we are using a four- vector notation where q = (q, u> n ) comprises both a wavevector and a Matsubara frequency, 
and J = T^2 n . V is the system volume, H denotes a magnetic field in z-direction. For notational simplicity we 
suppress the replica labels that are needed to carry out the disorder average [9. x^ an d X^ are disorder averaged 
correlation functions of a reference system of disordered electrons that has T t = jfl| . More generally, the coefficient 
of the term of order M" in the LGW functional is given in terms of an n-point spin-density correlation function 
in the reference system, and in writing Eq. (|^) we have left out terms which can be shown to be irrelevant for the 
critical behavior ||. We will come back to these omitted terms in the end. The particle-hole spin-triplet interaction 
r t is missing in the bare reference system, but it is generated by perturbation theory already at the one-loop level. 
The reference system then has all of the characteristics of the full system, except that it must not undergo a phase 
transition lest the separation of modes that is implicit in our singling out the spin-triplet channel for the H-S decoupling 
procedure break down. 

X^ (q) in Eq. (||) is the spin susceptibility of the reference system. At small Matsubara frequency and wavenumber 
it has the characteristic diffusive form [H, 
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X ( 2 )(q, W „) = X o(q)[l-K|/£q 2 ] 



(6) 



where D is the spin diffusion constant of the reference system, and Xo(q) is its static spin susceptibility. At T = 0, 
the latter can be shown to be a nonanalytic function of q of the form xo(q - * 0) ~ const — |q| d ~ 2 — q 2 , where we 
have omitted irrelevant prefactors (To). Choosing suitable units, the Gaussian part of <J> can then be written, 

$( 2 ) [M] = [ M(q) [t + c n \q\ d - 2 + c a q 2 + c d \io n \/q 2 } M(-q) , (7) 

Jq 

where to = 1 — I\x (<1 ~ > 0,w n = 0) is the bare distance from the critical point, and c„, c a and Cd are constants. 

The coefficient x^ m Eq. (JsJ) , contrary to the case of a usual LGW functional, is in general not finite at zero 
frequencies and wavenumbers. x^ is a nonlinear susceptibility in the reference ensemble. Such a quantity one 
expects to be divergent, given the singularity in the linear susceptibility Xo(q)- Indeed, standard perturbation theory 
shows H that it is given schematically by 

X (4) ~ const + v f[k 2 + \u n \] ~ 4 ~ u + v\p\ d - 6 . (8) 

Jk 

Here u and v are finite numbers, and we have cut off the singularity by means of a wavenumber |p|. The physical 
interpretation of this cutoff will be dealt with below. More generally, the coefficient of |M| 2 ™ in <i> behaves like 
t^-( 2 ™) ^ |p|d+2-4n^ ^his i m plies that $ contains a nonanalyticity which in our expansion takes the form of a power 
series in |M| 2 /|p| 4 . The nonanalytic structure of the correlation functions x*' 2 ™' 1 is a consequence of the diffusive 
spectator modes which we have integrated out. Therefore the M-ficld theory, rather than having a simple LGW form, 
is strongly nonlocal, which leads to unusual scaling behavior. 

The functional $ can be analyzed by using standard techniques ||. We look for a FP where Cd and either c n 
(for 2 < d < 4), or c a (for d > 4) are not renormalized. This fixes the critical exponents r\ and z. Power counting 
then reveals that both u and v are irrelevant. Several features of the critical behavior follow immediately from this 
observation. The critical exponents r\ and z are fixed by the choice of our FP, and v and 7 can be obtained by 
considering the q-dependence of the Gaussian vertex, Eq. (Q). We can thus read Eq. (Q) and the value of the exponent 
7 directly off Eq. (Q). 

We determine the equation of state by taking the term of order |M| 4 in $ into account. To this end, we need to 
identify the scaling behavior of the cutoff wavenumber p. Since p cuts off a divergence stemming from an integration 
over products of diffusion poles that is rendered finite by spin-flip processes, in the presence of an external magnetic 
field B one would identify p 2 with B. Since a nonzero magnetization acts physically like a magnetic field, it is therefore 
plausible that in the present context, and for scaling purposes, p can be replaced by to 1 / 2 . We have confirmed this 
conclusion by means of a more technical analysis The effective equation of state can therefore be written in the 
form of Eq. ( paf ) . 

These results completely specify the critical behavior at T — 0. Their most interesting aspect is the nontrivial 
exponent values found for 2 < d < 4, which can nevertheless be determined exactly. The reason for this is the |q| d_2 - 
term in the Gaussian action, Eq. (Q). It reflects the fact that in an interacting electron system with quenched disorder, 
static correlations between spins do not fall off exponentially with distance, but only algebraically like r -( 2d - 2 ). This 
leads to a long-range interaction that falls off like l/r 2d ~ 2 , see Eq. ( f7|). The critical behavior of classical Heisenberg 
magnets with such a long-range interaction has been studied before j^l) . 

We now turn to the T-dependences of the specific heat, cy, and the magnetization. Let us first consider the former. 
We expand the free energy functional, Eq. (|B|), about the expectation value, to, of M to second order, and then 
perform the Gaussian integral to obtain the partition function. The free energy is obtained as the sum of a mean-field 
contribution given by $[m], and a fluctuation contribution given by the Gaussian integral. The latter yields the 
leading nonanalytic term in the free energy. We find Q that effectively H and T have the same scale dimension, viz. 
d(= z), and that at t = there is a logarithmic T-dependence of 7^ for all 2 < d < 4. These are the logarithmic 
corrections to scaling mentioned above. If we put the t-dependence back in, we obtain that the scale dependence of 
7v is given by Eq. (^). 

For the magnetization, the situation is different compared to the specific heat, in that the leading T-dependence is 
given by the mean-field or saddle-point contribution to the free energy. We calculate the temperature corrections to 
the equation of state, Eq. (pa|), and find that for to >> T (in suitable units), m d ' 2 in Eq. ( pa| ) will be replaced by 
TO d//2 [l + const x T/to + ...], while for to << T, t is replaced by t + T 1 / 2 ". The effective scale dimension of T in m 
is therefore 2 (not z), and we obtain for m the homogeneity law given by Eq. (||). The appropriate interpretation of 
the relevant operator T in these equations is that it reflects the crossover from the quantum FP to the classical FP 
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rather than dynamical scaling. Accordingly, we have written the T-dependence in Eq. (g) in terms of the crossover 
exponent <p. 

We finally discuss some further aspects of our results, and their relation to other work. The phase transition 
described above is very similar to the one discussed in Ref. [ fl2| . In particular, the critical exponents on the disordered 
side and in the quantum critical region are the same, and Ref. pl| also found logarithmic corrections to scaling. We 
therefore believe that the quantum phase transition found in Ref. JL2| |, whose physical nature was unclear, was the 
ferromagnetic transition. Also note that throughout we have studied what amounts to the average effects of disorder; 
the coefficients in our LGW functional are disorder averaged correlation functions. In addition, fluctuation effects 
appear, e.g. in the form of a random mass term at 0(M 4 ). The coefficient of this term is the disorder average of the 
square of the unaveraged correlation function \^ ■ This term is irrelevant by power counting with respect to our FP 
for all d except d = 4, where it is marginal, and all other fluctuation terms are also irrelevant [^J. 

Sachdev [ p"3[ has performed a very general scaling analysis of quantum phase transitions with conserved order 
parameters. Some of his results do not directly apply to our problem because the coefficient v in Eq. (||) acts 
as a dangerous irrelevant variable (DIV). For instance, it follows from our Eq. (||) that the Wilson ratio, W = 
(m/ H) /(cy /T), diverges at criticality rather than being a universal number. The reason for this breakdown of 
general scaling is the flexibility displayed by DIV in modifying the scaling behavior of various quantities: Not only 
can a particular irrelevant variable be dangerous for some observables but not for others, it can also affect different 
arguments of a given scaling function in different ways. Here, the DIV v affects the magnetization differently than 
the specific heat coefficient, and in the latter it changes the effective scale dimension of H (from (3d — 2)/2 to d in 
2 < d < 4), but not that of T. 

We gratefully acknowledge discussions with Subir Sachdev and W.W. Warren. This work was supported by the 
NSF under grant numbers DMR-92-09879, DMR-92-17496, and DMR-95-10185. 
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